Slow holes in the triangular Ising antiferromagnet 
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We consider the problem of the doped Ising antiferromagnet on the frustrated triangular lattice 
in the limit where the hole kinetic energy is much smaller than the Ising exchange. For a single 
hole we prove a "frustrated Nagaoka theorem" showing that the ground state is magnetized and 
breaks translational symmetry, in contrast to the parent insulating state that is unmagnetized and 
spatially homogenous. The extension of this physics to finite dopings depends on the strength of a 
density-density coupling that is inevitably present — we find either phase separation of the holes, or a 
superconducting state that is also magnetized and breaks translational symmetry in a feat of spatial 
self-organization. Finally, we derive an effective interaction between dilute holes at temperatures 
in excess of the hopping and find an oscillatory, long-ranged form reflective of the correlations 
in the underlying classical magnet which presages the breaking of translational symmetry at zero 
temperature. 



PACS numbers: 75.10.Jm, 74.20.Mn, 05.50.+q, 71.27.H-e 



I. INTRODUCTION 

The problem of mobile holes in a magnetic background 
is an especially interesting example of charge motion in 
a strongly correlated electron system. In the cases that 
we have in mind, the magnetism involves local moments 
which reflect an underlying strong electron-electron in- 
teraction that can produce new and interesting effects 
when doping creates charge carriers. Much of the re- 
cent work on this problem has been inspired by the 
cuprate superconductors, whose physics is widely be- 
lieved to be centrally connected to their genesis as doped 
Mott insulators. El In this context, the central observa- 
tion is that hole motion in an ordered antiferromag- 
netic background, is frustrated p leading to the expecta- 
tion that dopinga will result in a rearrangement of the 
magnetic backgound in a manner conducive to reliev- 
ing this frustration. Suggestions for what this_ entails 
include RViB (resonating valence bond) theoryja stripe 
formation^ and, in models with purely short-ranged in- 
teractions, phase separation.!! A second theme in this set- 
ting is the large spin degeneracy of the extreme Mott in- 
sulator (e.g. the U = oo Hubbard model) at half filling 
and the role of doping in lifting it. Perhaps the most eel-, 
ebrated result along these lines is the Nagaoka theorem,Q 
which established that a single hole would completely po- 
larize the spin background. While the situation at finite 
dopings remains unsettled, the Nagaoka result does show 
that doping can lift the degeneracy in striking ways. 

In this paper we report some results on dilute holes in- 
troduced into a frustrated magnetic system — the particu- 
lar system studied is the canonical example of this class, 
the triangular lattice Ising Antiferromagnet first studied 
by Wannier and Houtappel.H This system realises aspects 
of both themes touched on above. It has local antiferro- 
magnetic order which leads to frustration of the hole mo- 
tion. In addition however, as this is an Ising system, the 
magnetic frustration leads to a finite zero-point entropy 
per site, S. This feature is reminiscent of the Mott in- 



sulator cited above, but we should note that in this case 
the degeneracy is generated as a cooperative effect: the 
Ising model on the triangular lattice has a large num- 
ber of ground states, NgSi because of the geometrically 
frustrated nature of its magnetic interactions. In each 
triangle, at least one pair of spins has to be aligned, and 
hence one bond frustrated. Any state with exactly one 
bond per triangle is hence a ground state, and Mgs is 
found to scale exponentially with the number of spins, 
N, with S = k B \n(Af GS )/N ~ 0.323fc B .B 

In the following, we present a set of results on the 
question of if, and how, this degeneracy is lifted upon di- 
lute doping. To make progress, we will assume the hole 
motion is slow — their kinetic energy being taken to be 
much smaller than the magnetic exchange. While Ising 
magnets exist, slow holes might be harder to find. Our 
interest in this limit then is that a) it poses the ques- 
tion of how doping interacts with a frustrated magnetic 
background most cleanly, which is of theoretical interest, 
and b) that the ordering patterns we find in this limit 
could well persist when more enterprising holes are con- 
sidered. While more work is needed to investigate the 
validity of (b) we should note concerning (a) that the re- 
sponse of highly frustrated magnets to perturbations is 
more generally interesting. On account of the large de- 
generacy, these systems are unstable in a large number 
of directions, promising a wide range of unexpected and 
unusual physical phenomena. Along these lines, in a re- 
cent study along with P. Chandra, we have explored the 
phase structure produced by switching on the auantum 
dynamics of a transverse magnetic field insteadJa 

We now turn to the Hamiltonian, H , we study. It is 

H = Hf + Hj + H v 
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where J z > is the strength of the antiferromagnetic 
Ising exchange, t is the hopping integral, and Cj CT (c| CT ) 
annihilates (creates) an electron at site i with spin a. 
Sf = |(cLcif — cj.Cjj) is the z-component of the spin- 

1/2 operator at site i, rii — ct is the total density 
operator and P is a projector that ensures that all sites 
are at most singly occupied. Finally, the sum (ij) runs 
over nearest neighbor bonds and the sum a over up and 
down spin orientations. As noted earlier, we study the 
range of parameters t <C J z , where the hopping acts as 
a perturbation to the exchange coupling. We consider 
strengths of the nearest neighbour density-density cou- 
pling of < rj < 1. In the derivation of the rotation- 
ally invariant t — J model for S = 1/2 from the Hubbard 
model, such a coupling is generated with rj = 1 . For Ising 
systems, possibly with spins greater than 1/2, there is no 
such precise fixing of its magnitude and for that reason, 
as well as for more general theoretical interest, we con- 
sider a range of values for it. 

In the rest of the paper, we first prove a "frustrated Na- 
gaoka theorem" , namely that the introduction of a single 
hole into the ground state manifold causes the system 
to order into a state that is spin polarized and breaks 
translational symmetry with a y3 x \/3-unit cell. In 
this state, the frustrated bonds form a hexagonal (honey- 
comb) lattice on which the hole hops. Next we consider 
generalizing this physics to a finite but dilute density 
of holes, which turns out to depend strongly on rj. For 
rj > 1/3, phase separation obtains. The case rj < 1/3 is 
much more interesting. Here, we find a spinless Fermi 
liquid state that lives on the "hexagonal backbone" of 
the single hole problem. Further, the surrounding sites 
mediate an attractive interaction that implies that the 
Fermi liquid is unstable to a superconducting state in a 
non-zero odd angular momentum channel. The resulting 
state then breaks three symmetries of the parent insulat- 
ing state at once: Ising, translation and the charge U(l). 
Finally, we turn to non-zero temperatures J z > T » t, 
and show using Grassman functional integral techniques, 
that there is an algebraically decaying, entropic interac- 
tion between two holes which is oscillatory as a function 
of distance and which indicates a tendency for the holes 
to segregate on the hexagonal backbone. 

II. A FRUSTRATED NAGAOKA THEOREM 

We show in this section that the ground state of our 
Hamiltonian (Eq. |l|) with a single hole present is macro- 
scopically ordered. We will begin by establishing this in 
the limit where we take J z — oo, i.e. in the degenerate 
perturbation theory problem among a set of single hole 
states that we identify below. We will comment on why 
the conclusion is unaltered when J z 3> t but not infinite, 
at the end of the proof of this first part. 

The proof will consist of first establishing a lower 
bound for the energy of a one-hole state on arbitrary 



finite lattices and then demonstrating that this can be 
saturated only on a subset of finite lattices and that on 
the latter there is, up to global symmetry operations, just 
one state that does so. For specificity, we will assume that 
our lattices come with periodic boundary conditions. 

We first identify the one-hole states that minimise the 
exchange energy, noting that the term H n in the Hamilto- 
nian (Eq. [l]) is the same for any one-hole state. Starting 
from any ground state, we can only remove an electron 
which experiences zero net exchange field, i.e., which has 
three frustrated and three satisfied bonds. There are ex- 
ponentially many such one-hole states. 

The hopping Hamiltonian connects those one-hole 
states which only differ by the exchange of a hole and 
a neighbouring spin. Fig. [lja shows a spin configuration 
which allows the hole to hop to three of its neighbouring 
sites. Up to symmetry operations, this is the only such 
(local) spin configuration. Since the hopping integral is 
t, this means that the energy which can be gained from 
H t is at most 3i. We note that this bound is quite gen- 
eral and holds for all finite lattices with arbitrary bound- 
ary conditions. Indeed, starting with any given one-hole 
state, we can think of the all the states reached by hops as 
inducing a graph, whose topology is in general far more 
complicated than that of the parent lattice. The state- 
ment is that regardless of the complexity, the graph has 
a maximum local coordination number of 3 and that the 
nearest neighbor hopping problem on it has at best an 
energy of — 3t. 

This energy gain can only be realised if any state 
reached after one hop again permits hoppping in three 
directions. Demanding that this be the case completely 
constrains the parent spin configuration to the one de- 
picted in Fig. |l|b; the other five states that satisfy this 
condition are obtained using translations or Ising rever- 
sal. In this configuration, the hole can occupy any site 
with a frustrated bond. With periodic boundary condi- 
tions, it is clear that not all finite lattices accomodate this 
pattern, but that there is a subset of arbitrarily large size 
which does — these are the lattices with linear sizes that 
are integer multiples of 3. In the following we construct 
our desired one- hole eigenstate \h) with energy — 3i on 
this subset. 

Define \G) = Y[i c l a (i) l^)> where a(i) is the spin at 
site i in the hexagonal state depicted in Fig. ^fa, and 
|0) is the empty lattice; the ordering of the operators is 
immaterial in what follows. Next, define the state with a 
hole on site n as \n) = {— l) s c„,j \G). Here, S(n) — (1) 
if the hole is created on sublattice A (B) of the bipartite 
hexagonal lattice defined by the frustrated bonds. The 
matrix elements (m\H t \n) vanish unless m and n refer 
to neighbouring sites on the hexagonal lattice, in which 
case 

(m\H t \n) = (m| - 4a c i<r W 
= -t (m\ c^ T c mT \n) 
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FIG. 1. (a) A spin configuration encircling a hole which 
allows the hole to hop onto the three neighbouring sites occu- 
pied by up spins. The spins are denoted by arrows, the hole 
is represented by 0. (b) The spin confiugration allowing the 
maximal coordination of three for a hole at any allowed site. 
Note that the sites of the up spins, on any of which a hole can 
reside, form a hexagonal lattice (the 'hexagonal backbone', 
fat lines) on which the hole can hop. 



can hop to a site where the net exchange field does not 
vanish and, depending on the spin configuration, either 
be forced to return to its starting site or hop to a different 
one. In terms of the induced graph mentioned above, the 
former corresponds to an on-site potential energy term, 
while the latter generates a new, or modifies the strength 
of an existing, edge (hopping amplitude). In the hexag- 
onal backbone states, only the first option is available, 
yielding an effective potential energy of —3t 2 /J z at each 
site. If we consider altering the state near a particu- 
lar site to generate more favourable potential or hopping 
matrix elements to 0([i/J z ] 2 ), this invariably requires 
sacrificing a hopping element of 0(t). We believe that 
this points to a certain rigidity of the hexagonal back- 
bone which will operate to exclude similar subtleties at 
higher orders as well. Needless to say, a proof of this 
assertion would be desirableU 



= -t(-l)(° +1 ) (G\cl n1 ci lCm1 c nl \G) 

= -t{G\c\ A c n ^ m] c m] \G) = -t . (2) 

Now, let \h) — Y^i={ 3 \ n )- For this state, we have 
H t \h) = —St \h) since H t connects, with amplitude —t, 
each single-hole component of the wavefunction to the 
three components with the hole on a neighbouring site. 
Note that this is, in essence, the ground state of a spinless 
hole hopping on a hexagonal lattice. 

Unlike the ensemble of undiluted ground states, which 
has neither long-range order nor a net magnetisation, 
| h) incorporates long-range spin order .El There is a three 
sublattice (v3 x y/S) ordering pattern with the spin on 
two sublattices pointing up and down on the third. The 
magnetisation induced by adding the single hole is there- 
fore AT/3. It is interesting to note that this ordering 
pattern would also have been generated by an explic- 
itly symmetry-breaking magnetic field pointing along the 
Ising axis and that the v5 x v3 structure is also picked 
out by a transverse field.L] 

Finally, we ask whether the conclusion at J z — oo is 
altered when J z is made finite. While we do not have 
a proof of this, we believe that this does not happen. 
Conceptually, it is useful to imagine generating an effec- 
tive Hamiltonian in the space of all one-hole states by 
including fluctuations outside this space perturbatively 
in powers of t/J z . Such an expansion will not invali- 
date the conclusion reached inside the hexagonal back- 
bone states — it will only modify the state \h) perturba- 
tively, preserving its symmetry characteristics at small 
t/J z . The real danger is that a state in a different sector 
that has an energy e = — 3t + 6(L) for a system of linear 
dimension L, when J z = oo, with limx^oo 9(L) = + , 
will acquire a perturbative gain of O ([t/ J z ] n ) for some 
n, in excess of the state \h), thereby invalidating our 
conclusion at sufficiently large L. 

To gain insight into why this is unlikely to happen, con- 
sider the energetics at O ([t/J*] 2 ) ■ At this order, the hole 



III. SMALL DENSITY OF HOLES AND 
LEADING t/.F CORRECTIONS 

We begin with the case r\ = 0. First, we need to find 
the allowed many hole states at low dilution (and when 
t/J z — > 0). As before, we need to remove those electrons 
experiencing zero net exchange field since this ensures 
that the leading term of the Hamiltonian, Hj, continues 
to be optimised. This implies that one cannot remove 
two neighouring electrons: this would leave a triangle 
with only one occupied site and hence no contribution 
to the magnetic energy, whereas the optimal states are 
those in which the three bonds of each triangle have total 
magnetic energy — J/4. 

We showed in the last section that the kinetic energy 
selects a subset of the one-hole states in which the hole is 
maximally delocalized. We believe that at sufficiently low 
dopings, this logic continues to operate and the low- lying 
eigenstates of H will be constructed out of the many-hole 
states obtained by removing electrons from the hexagonal 
backbone of the state depicted in Fig |l|b. (However, we 
are no longer able to prove this.) This state continues to 
be favourable as it allows the holes to be maximally mo- 
bile by having coordination three everywhere. The only 
restriction being placed on the holes is not to sit imme- 
diately next to each other, which is a demand present for 
all possible background states and one that should not 
be too onerous at low dopings. 

With these restrictions, the problem reduces to that 
of spinless fermions hopping on a hexagonal lattice with 
amplitude t and with a nearest neighbor hard core re- 
pulsion. At low densities, this is a-svell studied problem 
and known to be a Fermi liquid.113 In this fashion, we 
find that the ground state of the dilute hole problem is 
magnetized, breaks translational symmetry and is a two 
dimensional Fermi liquid. 

Now this is not quite true, for the Fermi liquid suf- 
fers from the well-known Kohn-Luttinger superconduct- 
ing instability!!^ and so strictly at T = the system will 
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break the charge U(l) symmetry as well — thereby break- 
ing three symmetries of the parent insulating state at 
once! While this happens as t/J z — > 0, we show next 
that for t/J z large but not infinite, the magnetic back- 
ground of the hexagonal conducting lattice mediates an 
attractive interaction that also drives a superconducting 
instability — at low densities this is a much larger effect 
than the purely "internal" Kohn-Luttinger effect. 

To this end we consider the effective Hamiltonian, 
H e ff, between holes on the hexagonal backbone to next- 
leading order in t/J z , induced by excursions to the sur- 
rounding sites in which the system leaves the ground 
state manifold of Hj. The matrix elements of this effec- 
tive Hamiltonian are determined by perturbation theory 
as follows. First, we define a fc-hole state with holes at 
sites mi, 7712, ... mu (with some ordering of the sites, and 
for concreteness with ro, < rrij for i < j) as 



TO1TO2 ••■ Tn r , 



II c *tI g > 



(3) 



1=1 



In the two-hole sector (k = 2) we thus obtain: 
{mim 2 \ H eff |mim 2 ) = ^ — So _ F o 



Pl,P2 



-; (4) 



m 1 ,r/!s Pi,P2 



here E° are the energies of the two-hole states for t = 
0. Note that the overall sign of the many-hole states is 
immaterial at this order as only the squares of matrix 
elements of Ht are used in determining H e f / . 

To this order, more general, non-diagonal matrix el- 
ements are absent since holes hopping of the backbone 
have to hop back to their original site in order to re- 
constitute a leading-order ground state. When the two 
holes are not on the same hexagon, H e ff always takes 
the same value which is just twice the contribution of an 
isolated hole and hence has the interpretation of an effec- 
tive one-body potential as discussed in the last section. 
Subtracting this leaves a term that depends on the joint 
presence of the two holes and can thus be viewed as an 
effective two-body potential, V(R), for the holes. This 
has the form, restoring 77 for future use: 



V(R = 1) =00 
V(R = VS)=- 

V{R = 2) = - 
V(R > 2) = 0, 



t 2 2(13 - 4?? -r/ 2 ) 
7 3 - 4?y + if 
t 2 2(1 + 7?) 
J 3-77 



where R is measured in lattice constants of the triangular 
lattice. 

This effective interaction is dominantly repulsive on 
account of the hard core piece considered earlier. Hence 
the additional, weak, attraction can only induce a super- 
conducting instability of the parent Fermi liquid and we 



do not need to worry about the possiblity of phase sepa- 
ration at 77 = 0. In the continuum, it is easy to see that 
the two body problem with a hard core and a weak at- 
tractive tail posesses attractive phase shifts in sufficiently 
high angular momentum channels — at low densities this 
is sufficient to lead to pairing. Ej While we have not car- 
ried out a detailed analysis on the hexagonal lattice, the 
general conclusion will hold. 

We end this section by considering the effect of switch- 
ing on an 77 > 0. For r\ < 1/3, the above arguments retain 
their validity, with the main consequence of rj being a 
strengthening of the attractive part of the effective inter- 
action. For 77 > 1/3, the allowed many-hole states change 
entirely in character. The density-density attraction H v 
then is strong enough to overcome the exchange term 
Hj and produce phase separation. This happens because 
phase-separation costs an energy of J(3y7 + 1)/4 per hole, 
whereas the optimal states described above achieve a cost 
of 3J77/2. 



IV. EFFECTIVE INTERACTION FOR 
CLASSICAL HOLES 

Thus far we have focussed on the physics precisely at 
T = 0. In this limit the degeneracy of the parent magnet 
is "selected away". At non-zero temperatures, but still 
below J z , the ground-state entropy of the parent magnet 
will again play a role. As a first step towards an under- 
standing of finite temperatures we consider t <C T <C J, 
equivalently, we discuss the behaviour of classical an- 
nealed holes on the triangular antiferromagnet.tJ Our 
central result in this part is an effective interaction, in- 
duced entropically, between two holes. We close with 
some remarks on the finite density problem. 



A. Two holes 

Two holes on the triangular lattice experience an en- 
tropic interaction because the number of ground-state 
spin configurations, Z(r), minimising Hj depends on 
their separation vector, r. From a knowledge of Z(r), 
upon subtracting the "creation entropy" of the two holes, 
ln.Z'(oo), we can obtain an effective interaction potential 
(3v(r), which vanishes as r — > 00: 



(3v(r) — — In Z(r) + In Z(oo) 

= - \n{Z(r)/Z) + ln(Z(oo)/Z); 



(5) 



here Z is the partition function of the undiluted system.t3 
Defining Z{oo)/Z = T 2 and Z(r)/Z = T 2 - ((r), we 
obtain 



I3v(r) = — In 



1 - 



CM 

y 2 



J~2 ' 



(6) 



the last step being justified at large \r\ by the smallness 
of £(r)/Y 2 = 0, as we will show below. 
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To determine ((r) and T, it is convenient to use a 
standard representation of the spin problem as a clas- 
sical dimer model on the dual hexagonal lattice. The 
presence of a hole is encoded by a certain dimer con- 
figuration around it. The latter is soluble by the Pfaf- 
fian techniques introduced by KasteleynEZl although we 
will find it convenient to use the language of Grassma- 
nian functional integrals introduced into such problems 
by SamuelJla Thus, the function £(r) is obtained by eval- 
uating a twelve-fermion correlation function. 

In detail, this calculation proceeds as follows. First, 
we map ground state spin configurations on the triangu- 
lar lattice to dimer configurations on its dual hexagonal 
lattice by drawing a dimer through each frustrated bond, 
placing its endpoints at the centres of the two triangles 
sharing the frustrated bond. Noting that each triangle 
has one and only one frustrated bond in the ground state, 
the dimers thus provide a hard-core covering of the lat- 
tice dual to the triangular lattice, namely the hexagonal 
lattice. Conversely, up to the global Ising symmetry, each 
dimer covering corresponds to a ground state of the Ising 
model. 

The statistical mechanics of classical hard core dimer 
coverings-en the hexagonal lattice has been studied by 
Kastelynta and some correlation functions have been 
computed by Yokoi, Nagle and Salinas.Ell 

We extend this method to include the presence of di- 
luted sites on the original triangular lattice. As outlined 
above, a hole can only occupy the site of a spin experienc- 
ing a net exchange field of zero. In dimer language, this 
means that only a spin at the centre of a hexagonal pla- 
quette occupied by exactly three dimers can be replaced 
by a hole, as depicted in Fig. |^. To determine Z(r), we 
arc thus only allowed to count those dimer configurations 
in which such plaquettes are encountered at the location 
of the holes. 

To calculate Z in the first instance, one defines Grass- 
man variables, ipi, on each site of the s sites on the hexag- 
onal lattice. One then writes 



Z = 



J fn^'J etPtMrti) = Pf ( A ) = V / deII, (7) 



where Pf(-A) is the Pfaffian of the (antisymmetric) matrix 
A. 

The problem of determining A for a general class of 
two-dimensional lattices has been solved by Kasteleyn. 
For the hexagonal lattice, which we are interested in, it 
turns out that one has to double the unit cell of the lattice 
to contain 4 spins, with the unit cell labelled as in Fig. [|, 
so that Aij(r) equals 



/ S r>0 - 6 r> £ 

— S r> o + 5r,-£ 
V S r ,-v 



S r _ v 

S rfi - 5 r £ 

S r fl + <5 r ,_£ 

Here, the vectors £ and v denote the translation vectors 
of the lattice with the doubled unit cell (Fig. 0). The 




FIG. 2. The hexagonal lattice with the spins of the unit 
cell labelled 1-4. The lattice translation vectors are denoted 
by £ and v for the x and y directions, respectively. The pla- 
quettc on the top right encloses a hole, marked by a cross. 
One of the two possible resulting dimer plaquette configura- 
tions is shown. 



holes are located in cells r±, site i, and ra, site j, with 
r = r-z — ri . 

This choice of A accomplishes the following. If the ex- 
ponential in the integrand for the partition function is 
expanded, the only terms that result in a nonzero contri- 
bution to the integral are those in which each ipi appears 
exactly once. Each such term represents a covering of 
dimers with sites i and j being connected by a dimer if 
Aij is one of the prefactors of the term. The signs of the 
entries of A are chosen such that every such term in fact 
integrates to 1 so that all dimer coverings are allocated 
equal weight. One thus obtains S/ks = In Z/N ~ 0.323. 

We first caclulate the probability that a spin can 
be replaced with a hole, which we require as this will 
turn out to be T, defined above. To ensure the pla- 
quette surrounding a spin (the sites of which we label 
1...6) is covered by three dimers, we insert a pref ac- 
tor AyiA^A^ilix^i^^^}^}^ into the integrand, which 
forces the three dimers into the required position. We 
thus need to calculate 



T 




I x 



Ai2A 3 4A 56 ip 1 ip 2 4>34>4i>5i>6 exp^Aij^j). (8) 



Since the Fermionic action is quadratic, this becomes 
by Wick's theorem the sum over all 6! contractions of 
pairs of ips: 



T = 



^ sign(P)Gi 1 i 2 Gi 3 i 4 Gi 5 i s , 



(9) 



where the sum runs over all permutations P = = 
1...6} of {1...6}. These Green functions have been worked 
out by Yokoi et altQ together with identities expressing 
all Gij in terms of G24, which for y > is given by 



r 77 / 3 

g m (x,v) = -(-ir 



cos(2a;0) 



|l + exp(2i<£)| 



2(»+l) 



(10) 
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with G24 given by the negative of this expression for y < 
0. 

In fact, the sum for T only contains 3! nonzero terms 
as Green functions arising from contractions of sites be- 
longing to the same sublattice of the hexagonal lattice 
vanish. We thus obtain 



2 

27 



3\/3 

"8^3" 



1 



2a/3 tt 



0.14501. 



(11) 



To confirm the calculation so far, we remark that Monte 
Carlo simulations we have carried out in a different 
context (frustrated transverse field Ising models)B give 
T = 0.147 ±0.003. 

The evaluation of Z(r) thus requires twelve-fermion 
correlation functions as two holes have to be introduced. 
There are two technical difficulties in this calculation. 
Firstly, the distance-dependence of the Green functions 
is not known in closed form, and secondly, we now require 
6! = 720 terms each containing at least six Green func- 
tions even after using the vanishing of same-sublattice 
Green functions. 

As the reader can see, this calculation is straightfor- 
ward in principle but somewhat involved in practice — it 
will turn out though that the final answer has a trans- 
parent rationale. Before proceeding further it is useful 
to discuss the same general phenonmenon in a toy model 
which is exactly solvable and where the combinatorics 
is more transparent. Consider a one-dimensional dimer 
model in which the dimers cover a two-leg ladder (see 
Figf|. The entropy per rung of this dimer model is the 
golden mean, G = (a/5 + l)/2 ~ 1.62. The dimer model 
on this ladder arises ia the study of a fully frustrated 
three-leg Ising ladder, til and removing a spin from this 
Ising model corresponds to forcing the square surround- 
ing this spin to be covered by two parallel dimers. 

If two holes are adjacent to one another (Fig. ||b, top), 
there are no rungs separating the plaquettes surround- 
ing them. Moving the hole on the right two units fur- 
ther away (Fig. |^b, bottom) costs an entropy of InG 2 ~ 
In 2.62 in the semi-infinite region outside to the right of 
the hole pair and only gains an entropy of In 2 as the 
dimer pair between the holes can only have two states. 
For even separations, it is thus advantageous for the holes 
to move close together. 

Separating the holes by one unit costs entropy InG 
outside while gaining nothing inside as a single dimer can 
only exist in the vertical state (Fig. ||b, middle). Moving 
the hole further away by two units, however, yields a 
total gain of In 3 — InG 2 > in entropy as the three 
dimers between the hole plaquettes can have three states. 
Therefore, for odd separations, the holes want to move 
further apart. In fact, the alternating two-hole potential 
for this model can be calculated both by Pfaffian and 
transfer matrix methods, and takes the form 
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FIG. 3. (a) The fully frustrated three-leg ladder. Dashed 
(solid) lines indicate antiferromagnetic (ferromagnetic) ex- 
change. The dimer state corresponding to all spins pointing 
up is shown, (b) The corresponding dimer model on the dual 
lattice. A hole, marked by a cross, fixes the dimer plaque- 
tte surrounding it. Odd hole separations are entropically un- 
favourable because dimers can only resonate in pairs (bottom 
graph) and not on their own (middle graph) . 



with T ladder = 1/ (G (2G - 1)). 

Returning to the triangular lattice problem, the large 
|r| asymptotics of the Green functions are calculated in 
Rcf. M. 



G M {x,y) 



V3y cosm^) + x sin(^p; 

TT X 2 + 3 TT y 2 



(13) 



As these Green functions decay algebraically with |r|, 
the terms can be arranged according to how many con- 
tractions between sites belonging to different holes they 
contain. The leading terms are the ones in which all 
contractions are between sites of one or the other hole, 
and these terms are clearly independent of r. In fact, 
it is easy to see that the contractions around a given 
hole are simply those which yielded T in the single-hole 
calculation, and thus the leading, r-independent term is 
Z(oo)/Z = T 2 , as promised above. 

The leading non-trivial r-dependence arises from terms 
which contain two Green functions connecting sites 
around the two holes, an odd number of connections be- 
ing impossible. There are 324 such terms. With the 
assistance of Mathematica, we were able to analyze their 
asymptotics and obtain the exact long distance form, 



_cos(47rx/3) 
x 2 + 3y 2 ' 



(14) 



where 



729 - 324 V3 tt - 324 tt 2 + 96 \/3 tt 3 + 64 tt 4 



288 tt 6 



0.025853. 



(15) 



Translating the co-ordinates back to the original trian- 
gular lattice, our final asymptotic expression for the in- 
teraction potential is 



(3v(r) 



3 cos(47ra;/3) 
T 2 x 2 + y 2 



(16) 



The physics of this expression is quite transparent. 
The wavevector q = 4^ x on the triangular lattice is 
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characteristic of the hexagonal frustrated bond pattern 
sketched in.Fig[j]b. It has long been recognized (already 
by Wanniem) that this pattern represents the dominant 
ordering tendency in the ground state average of the clas- 
sical problem. In a moxe modern formulation in terms of 
height representations^ the pattern is the reference flat 
surface about which the system exhibits critical fluctua- 
tions. In accordance with this we find from Eq. [l^ that 
the holes attract when they share a hexagonal backbone 
of frustrated bonds, and repel when they do not. The 
algebraic decay could have been anticipated in the dimer 
language from our observation that the leading piece at 
large distance comes from a sum of dimer-dimer corre- 
lation functions each of which individually decays as the 
inverse-square of the distance. However, and this is what 
motivated the computation, it is not evident in advance 
that the various dimer correlators (which enter with dif- 
ferent signs) do not sum to zero and hence the acutal 
result does not just yield a non-trivial constant but also 
assures us that it is not zero. 



B. Finite density of classical holes: an 
incompressible state 

As the two hole interaction derived previously has con- 
siderable spatial structure and a long range, the reader 
may wonder whether it leads to a condensation of a finite 
density of holes on the hexagonal backbone. This does 
not happen. To see this note that f3v(r) is independent 
of T, which is to say that there are flucutations even as 
T — > — as indeed there are on account of the degeneracy. 
In the dilute limit, in which the two hole calculation is 
reliable, the effective interaction strength is proportional 
to the inverse square of the hole separation, i.e. to the 
(low) density and hence the hole gas remains in a liquid 
phase. Nevertheless, the long range of the interactions 
indicates that it exhibits substantial correlations coming 
from the entropic potential, and — on account of the form 
of the interaction — we suspect that such correlations will 
be algebraic as well. As we have already argued that the 
zero temperature state involves hole condensation on the 
hexagonal backbone driven by the hole kinetic energy, it 
appears that this phenomenon arises for both energetic 
and entropic reasons, and is thus quite robust so that it 
will characterize the entire temperature range T -C J z . 
(In all of this we are again referring to 77 < 1/3, else 
the energetic attraction ignored in the entropic deriva- 
tion will again dominate and drive phase separation.) 

Finally we report a curiosity regarding classical holes 
for rj < 1/3. As the density of holes is increased, it is 
not clear whether the interaction finally becomes strong 
enough to produce phase separation into regions either 
fully occupied or maximally diluted (but respecting the 
hard core repulsion of the holes). However, even if it does 
not, there exists a density of holes at which long-range 
order is present. 



The closest packing compatible with the ground-state 
constraint is obtained by placing the holes such that they 
cover one of two sublattices of the hexagaonal backbone 
of up spins depicted in Fig. 0b. At a hole density of 1/3, 
this sublattice is completely occupied by holes, at zero 
energy cost at O(J) compared to the undiluted system. 

Adding further holes above this density costs a mag- 
netic energy of order J, and thus the chemical potential 
experiences a discontinuity at a hole density of 1/3. The 
state at this density, which has perfect unfrustrated an- 
tiferromagnetic order on the bipartite hexagonal lattice 
occupied by the remaining spins, is thus an incompress- 
ible state. 

As the density of holes is increased further, phase sep- 
aration occurs for certain. Completely depleted regions 
are created by removing one hexagonal ring after another 
as the system optimises its energy by removing those 
spins with the smallest number of neighbours. This is 
the same mechanism that would operate in any unfrus- 
trated magnet. As mentioned above, for 77 > 1/3, phase 
separation occurs for all hole densities. 

V. SUMMARY 

We have considered the problem of dilute, slow holes in 
the triangular Ising antiferromagnet. We find that dop- 
ing has two very different outcomes depending on the 
dimensionless strength 77 of the density-density interac- 
tion. For 77 < 1/3 we find quite generally that at all 
temperatures T <C J z the holes tend to condense on 
a hexagonal backbone of frustrated bonds and that at 
T = they form a superconducting state that coexists 
with magnetic order and the breaking of translational 
symmetry in a "magnetic supersolid" . As a special case 
of this more general phenomenon, we were able to prove 
a "Nagaoka" theorem for a single hole at J z —00. For 
7/ > 1/3 in our short ranged model, we find phase sepa- 
ration. Whether the inclusion of the long range piece of 
the Coulomb interaction will change that outcome and 
stablize the magnetic supersolid and how much of this 
structure will persist to larger values oft/J z remain top- 
ics for future work. 
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